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©arge Order Behaviour of 2D Gravity Coupled to d < 1 Matter 
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We discuss the large order behaviour and Borel summability of the topological ex- 
pansion of models of 2D gravity coupled to general (p, q) conformal matter. In a previous 
work it was proven that at large order k the string susceptibility had a generic a k Y{2k— -|) 
behaviour. Moreover the constant a, relevant for the problem of Borel summability, was 
determined for all one-matrix models. We here obtain a set of equations for this constant 
in the general (p, q) model. String equations can be derived from the construction of two 
differential operators P, Q satisfying canonical commutation relations [P, Q] = 1. We show 
that the equation for a is determined by the form of the operators P, Q in the spherical or 
semiclassical limits. The results for the general one-matrix models are then easily recov- 
ered. Moreover, since for the (p, q) string models such p = (2m + l)q ± 1 the semiclassical 
forms of P, Q are explicitly known, the large order behaviour is completely determined. 
This class contains all unitary (q+ 1, q) models for which the answer is specially simple. As 
expected we find that the topological expansion for unitary models is not Borel summable. 
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1. Introduction 



We report here new results concerning the large order behaviour of the perturbation 
series of models [[]]-[§ of 2D gravity coupled to D < 1 matter. Our motivation is to 
gather some information about non-perturbative features of quantum gravity and string 
theory studying the asymptotic behaviour of perturbation series. In particular we want 
to understand whether perturbation theory indeed provides a proper definition of the 
physical theory of interest, in more technical terms whether the perturbation series is 
Borel summable. 

We first recall that the coupled differential equations for the partition function in the 
formulation of 2D quantum gravity coupled to arbitrary (p, q) minimal conformal matter 
can be derived from canonical commutation relations [P, Q] = 1 |] where P, Q are two 
differential operators of degree p and q respectively: 

p = dP - \ dP-2l >' Q = dq - 1 E w*)' d<? ~ 2 ^ ( L1 ) 

i=l i=l 

and u(x) = U\{x)/p = vi(x)/q is the specific heat or string susceptibility. Note that our 
normalization of u(x) differs by a factor 2 from the most commonly used in this problem 
(this normalization corresponds in the one-matrix case to consider potentials which are 
not even). In this way the double pole of smallest residue of u(x) has residue 1. Since the 
partition function F is given by F"(x) = —u(x), e F has then simple zeros. 

When one of the operators is given it can be shown that the other operator can be 
taken of the form: 

p = Qi q 

where the subscript + means that P is the sum of the terms of non-negative power in the 
formal expansion of Q p / q for d "large". In § (see also §), it was shown that the coupled 
differential equations also follow from an action principle. The basic action for a critical 
(p, q) model takes the general form 

S = J dx (Res Q p / q+1 + im) , (1.2) 

where Res denotes the residue (coefficient of <i _1 ) of its fractional powers. 

In the simple one-matrix case the "string equation" for the specific heat u{x) [0-0 
reduces to: 

(I + ±)Ri[u] = x , (1.3) 
1 



where the R^s are the usual KdV potentials ]7J. Due to the elementary properties of the 
i?z's, the above equation follows as the variational derivative with respect to u of the action 

S = J dx (Ri+i[u] +xu) . (1.4) 

In the following sections, we shall combine these properties with a direct analysis of the 
differential equations satisfied by the partition functions of the d < 1 models to determine 
the large order behaviour of the topological expansion of their solutions. Previous work || 
has allowed to determine that the topological expansion of the specific heat had the general 
property of behaving like a k T(2k — l/2) for k, the order in the topological expansion, large. 
The constant a was determined as the solution of an explicit algebraic equation for the 
one-matrix model (q = 2) and in two examples the critical and tricritical Ising model ((3, 4) 
and (4, 5) models). The importance of an explicit determination of a relies on the following 
property: If a is real and positive the perturbative expansion is not Borel summable and 
does not determine a unique function. Moreover there are good reasons to expect the 
corresponding model to be actually unstable. Such a result was obtained for half of the 
one-matrix models (this includes pure gravity), and is expected for all unitary (q + l,q) 
models. This latter property is derived here and the more general models p = (2m + l)g±l 
are explicitly discussed* . 



2. Large order behaviour of pure gravity 

We first recall the derivation of the large behaviour of pure gravity, because it illus- 
trates several features of the general analysis. 

For pure gravity, the differential equation satisfied by u(x) is 

u 2 (x) - -u"(x) =x. (2.1) 
6 

If u(x) has an asymptotic expansion for x large, it satisfies u(x) = ±y/x + O (x~ 2 ). The 
solution that corresponds to pure gravity has a x large expansion of the form 

u(x) =x 1 / 2 (l-Y,u k x- 5k ' 2 ) , (2.2) 

k=i 

* For a recent treatment of some standard features of divergent series, Borel summability, and 
summation methods, with physical applications, see, for example, pp. 840-842 of ||. 
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where the Uk are all positive. 

To determine the large order behaviour of the expansion we first analyze the stability 
properties of the solution for x large. Let us set u(x) i— > u(x)[l + e(x)] in eq. (|2.1] ) and 
write the equation obtained by expressing that the term linear in e vanishes: 

12u- — )e-2-e' -e" = 0. (2.3) 
u ) u 



One verifies that at leading order for x large only the leading order in eq. ( |2.2|) is needed 
and u" /u is negligible. Eq. (|2.3|) can then easily be solved by the WKB method. We set 



e'/e = r^i + bu'/u + O (u ,2 /u 5/i 



and find r = 2\/3 and b = —5/4. Replacing u by its asymptotic form u ~ x 1 / 2 and 
integrating we obtain: 

8^^.5/4 

e{x) oc x" 5/8 e 5 . (2.4) 



To leading order, the function e is also proportional to the difference between any Borel 
sum of the series and the exact non-perturbative solution of the differential equation (up 
to even smaller exponential corrections corresponding to multi-instanton like effects). In 
terms of the expansion parameter (string loop coupling) k 2 = x -5 / 2 , e reads 



8 (v/3> 



(x(k)) oc K 1/2 e 5 . (2.5) 



The above solution is valid for x large, i.e. k small. The large order behaviour in ( |2.2| ) is 
then given by 



d/; , , / 5 x 2 ' 



(The constant of proportionality in the above cannot be determined by this method.) The 
asymptotic T(2k — -|) behaviour is a slight refinement of the (2/c)! behaviour determined 

in HBjrg. 

The reality of r 2 has implied that all terms at large order have the same sign. This 
induces a singularity on the real positive axis in the Borel plane, obstruction to Borel 
summability. 

In it is confirmed that the exponential in ( |2.4|) coincides with the action for a 
single eigenvalue climbing to the top of the barrier in the matrix model potential, allowing 
us to interpret the exponential piece of the solution to ( pTT| ) as an instanton effect. 



3. The general string equations 



3.1. The general one-matrix problem 

We now consider the string equation ( |1.3|) , Ri[u] oc x. Substituting as before u(x) h- > 
u(x)(l + €(x)) we get a linear equation for e. At leading order for x large we expect the 
equation to be again solved by the WKB ansatz e'/e = rit 1 / 2 . It is then easy to verify that 
to obtain the leading large order behaviour of perturbation theory, it is only necessary to 
know the terms in Ri[u] that contain at most one derivative of u factor. The next leading 
contribution is given by terms such as ui-V 21 - 2 *" 1 )^, i.e. with a single factor of u' as 
well. 

At leading order only the terms in which the derivatives act on e are relevant and thus 
e satisfies an equation of the form 

I 

= J2 A ij uj ~ 1 z (2l ~ 2j) ■ (3- 1 ) 

3=1 
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The WKB ansatz leads to an (I — l) st order equation for the constant r 

= Mr) = J2 A » r2l ~ 2j - ( 3 - 2 ) 

3 = 1 

From the properties of the R^s one can derive an explicit expression for the polynomials 

where the subscript + again means the polynomial part of the large r expansion. The 
function (z 2 — l)^ 1 ^ 2 is also proportional to C^i+iiz) where C% 1+1 is a Gegenbauer poly- 
nomial defined by analytic continuation in v |§ . Note that the number of zeros is exactly 
the same as the number of operators in a (p = 21 — 1, 2) minimal conformal model |[l2[ . 
This is a property we shall meet again in the general case. Actually in the one-matrix case 
there is a natural explanation for it. The steepest descent analysis shows that the number 
of different instantons is related to the degree of the minimal potential corresponding to a 
critical point. This degree in turn is also related to the number of relevant perturbations 



We thank F. David for this remark. 
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For / even, eq. ( |3.2|) is an odd-order equation that will have at least one real solution 
for r 2 , positive as is obvious from the integral representation ( |3.3j ). The series therefore 
cannot be Borel summable. 

For I odd, on the other hand, the equation ( |3.2| ) for r 2 has no real solutions and there- 
fore we expect the solution of the differential equation to be determined by the perturbative 
expansion. 

Actually there exists a direct correspondence between the property of Borel summa- 



bility and the existence of the original integral. It has been noted |T3[ that according to 
whether I is odd or even, the original minimal matrix integral is well-defined or not because 
the integrand goes to zero in the first case while in the latter case it blows up for M large. 
A direct calculation, using steepest descent, of the instanton action || confirms that when 
the potential is unbounded from below the instanton action is indeed real and the series 
therefore non-Borel summable. 

The subleading terms in Ri[u] are immediately deduced from the leading terms by 
noting that since Ri[u] is derived from an action, eq. (|1.4|), the operator acting on e is 
hermitian. Therefore the operator w- 7_1 d 2Z_2: ' should be replaced by the symmetrized 
form d correcting ( |3.1| ) to 

I 

= J2 A ij{^~^ {2l ~ 2j) + H 21 - 2J')(J + ly-VeO"- 2 '- 1 )) . (3.4) 

3=1 

To characterize more precisely the large order behaviour, to next order we set e'/e = 
ru 1 / 2 + bv! /u, from which it follows, to the same order, that 

l—= r k u k/2 +r k-l u (k-3)/2 u , k ( b+ l(fc_l)) . (3.5) 



Substituting into ( |3.4|) , we find b = —{21 + l)/4, independent of r. Then 

21 (21+D/21 



e(x) oc x-( 2l+1 ^e 21+1 . (3.6) 

generalizing ( |2.4| ). In terms of the expansion parameter k = a-^+D/ 2 ', we find 

2 k 

«* A./^ e W"» K (^) r(2fc -^»' (3 - 7) 

The T(2k - \) factor in ( |2.6| ) is general, and is related to the property that the original 
equations descended from an action principle. 
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3.2. The general (p, q) model 

In the case of the general (p,q) model (eq. (|1.2| )) there results a system of coupled 
linear differential equations for the variations e Ui (x) associated with the functions Ui(x). 
At leading order we set e' u ./e Ui = to 1 / 2 . We obtain, taking into account the leading 
relations between the Ui, a linear system for each of the e Ui 's multiplied by a power of u 
determined by the grading. Imposing again the vanishing of the determinant of the linear 
system gives an equation for the coefficient r (and to leading order all functions e Ui are 
thus proportional up to a power of u determined by the grading) . 

To determine more precisely the behaviour of e u = e Ul we have to consider subleading 
terms. As in the one-matrix case they can be determined by a hermiticity argument. 
Since the equations for Ui derive from an action (|1.2| ), the linear equations for e Ui define a 
hermitian operator. This property leads to a universal T(k — |) behaviour for all the (p, q) 
models. 



We recall finally that the (2k)l large order behaviour is also the generic behaviour [|14 
for d = 1 models f|l4l 



4. The p = (2m + l)q ± 1 models in the spherical limit 

In the analysis of the large behaviour the knowledge of the solutions of the string 
equation in the spherical limit was required. Actually we shall prove in next section that 
the knowledge of the differential operators P, Q in the same limit is sufficient. For q = 2 
and q = 3 the form of the operator Q and P is known. In the general case q > 4 the explicit 
functional form of the operator Q in the spherical limit depends on the (p, q) models. In 
particular in the spherical or semiclassical limit the operator Q takes the form 

Q(d,u) = y £ i q i u i (x)&*-™ (4.1) 

i=0 

but the coefficients qi in equation ( |4.1| ) are in general p-dependent. Note that, in this limit, 
the order between the operators u(x) and d is irrelevant. 

However in [|16| it has been shown, using the string actions (|1 . 2| ) [0 in the semiclassical 
limit, that when p = (2m + l)q ± 1, the semiclassical form of the operator Q is m- 
independent and P and Q can be determined explicitly. This property can be recovered 
by a direct method. If we also set: 

P(d,u) = Y,P*u l (x) ( F- 2 \ (4.2) 

i=0 
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we obtain the semi classical limit of the equation [P, Q] = 1: 



,(dPdQ dQdP\ 



v <9d dit <9d <9it 
We now use the homogeneity property of P, Q: 

P(d, u) = u p / 2 P (du' 1 / 2 , l) , Q(d, m) = W 2 Q (dw" 1 / 2 , l) . 

From now on we call -P(^), Q(^) the two polynomials P(z = dw. -1 / 2 , 1), Q(z = dw -1 / 2 , 1). 
They thus satisfy the differential equation: 

qP'(z)Q(z) - P P(z)Q'(z) = 2pq , (4.4) 

while as expected the equation for u(x) yields 

u (p+q-l)/2 K x _ 

When one of the polynomials 

is known the other is obtained by integrating the equation. In the special case p = 
(2m + l)q ± 1 the polynomials Q(z) are Tchebychev's polynomials. Setting z = 2 cos 9 one 
finds that Q(z) = 2 cos q6 and 

P = 2pQ p / q (z) [ Q^-P/i^dt oc 2 V MA cos(p - 2lq)6 , 
Jo l=0 V I J 



satisfy the equation. 



5. Instantons: A more direct method 

5.1. Instantons in the one-matrix model revisited 

Before discussing the general unitary models let us return to the one-matrix model 
for which the result is exactly known. From the analysis of the corresponding non-linear 
differential equations we have learned that if we call e the variation of the specific heat 
u(x) then it has for x large the asymptotic form: 

e'/e~r\/u, (5.1) 

where r is constant which is determined by an algebraic equation. Since the variation of 
u can be neglected at leading order, we can rescale d, i.e. set m to 1. Equation ( |5.1| ) can 
then be written as a commutation relation 

de = e(d + r) /(d)e = e/(d + r). (5.2) 
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Then the operators P, Q are simply 



Q = d 2 -2, P = P 2m (d) = (d 2 -2)^. 

The equation for e is obtained by expanding at first order in e the commutation relation 
[P,Q] = 1. Setting: 

5P = { e , J R(d)} = ^i? fc {e,d 2 '- 1 - 2fc }, 

A;=0 

we find: 

[{e,P(d)},d 2 -2] + [P,-2e]=0. 
Using the commutation relation (|5.2| ) to commute e to the left we find the equation: 

-(2rd + r 2 ) (P(d)+P(d + r)) -2(P 2Z+1 (d + r)-P 2Z+1 (d)) =0. 

The first term vanishes for d = — r/2, which must thus be a zero of the second term. 
Taking into account the parity of P 2 z+i we obtain 

P 2 ; +1 (r/2)=0, (5.3) 

in agreement with the direct calculation. The polynomial P(d) is then determined by 
division. 

5.2. General (p,q) problem 

In the general (p, q) case, in the same classical limit, and after the same rescaling we 
have: 

Q = Q(d), P = P(d) = Q^(d), 5Q = {S(d) 1 e} 1 5P = {P(d),e}, 

where P, Q are polynomials of degrees p, q respectively, and P, S of degrees p — 2, q — 2 
and same parity as P, Q. 

The equation for e then leads to 

[P, SQ] + [SP, Q} = 
(P(d + r)-P(d))(S(d) + S(d + r)) 
- (Q(d + r) - Q(d)) (P(d) + P(d + r)) = . 

8 



The polynomials P(d+r) — -P(d) has a degree p—1 in d, while R has only a degree p — 2. An 
equivalent property is true for Q, S. The polynomials P(d + r) — P(d) and Q(d + r) — Q(d) 
must thus have at least one common root. Note that the first polynomial has p — 1 roots 
and the second q — 1. Moreover this roots are symmetric in the exchange d i— > — r — d. 
Therefore expressing the existence of a common root leads to (p — l)(q — 1) values of r, up 
to the symmetry. Note that the number of zeros is again exactly the same as the number 
of relevant operators in a (p, q) minimal conformal model [[12] of gravitationally dressed 
weights 

p + q-\pn-qm\ 

dm,n = ; ; , (5.4) 

p + q — 1 

with 1 < n < q — 1, 1 < m < p — 1 with the symmetry d m ^ n = dq-n^-m. The explanation 
of this relation is probably again that the number of different instanton actions is related to 
the degree of the minimal potential needed to generate a critical point in the multi-matrix 
model, and thus to the number of different relevant operators. Also we note that we are 
studying a general deformation of a critical solution and therefore the appearance in some 
form of the relevant operators should be expected. 

This condition determines the possible values of r when the polynomials P and Q, i.e. 
the differential operators are known in the classical limit. Examples are provided by the 
models p = (2m + l)q ± 1 where integral representations for these polynomials have been 
found. The simplest examples are provided by the q + 1, q models, i.e. the unitary models 
which we examine below. 

Let us finally verify that we can then indeed find the polynomials P, S. We call a 
the common root and assume first that a ^ — r/2. Then the parity properties imply that 
— r — a is also a common root. Setting then 

(P(d + r) - P(d)) = (d - a)(d + r + a)P(d), 
(Q(d + r) - Q(d)) = (d - a)(d + r + a)Q(d), 

we find that R and S are solutions of: 

R{d) + R{d + r) = (d + r/2)P(d), S(d) + S(d + r) = (d + r/2)Q(d). 

Note that these equations satisfy both the degree and parity requirement. 
If a = —r/2 the situation is even simpler 

P(d) + R(d + r) = (P(d + r) - P(d))/(d + r/2) 
S(d) + S(d + r) = (Q(d + r) - Q(d))/(d + r/2). 



9 



5.3. The unitary models 

We have shown that the differential operators P, Q may be written in the classical 
limit as: 

P = 2T p (d/2), Q = 2T g (d/2), 
where T p is the p-th Tchebychev's polynomial: 

T p (cos</?) = cospip. 

As explained above, taking into account the degrees of the polynomials R and 5, we 
conclude that the polynomials T q ((r + d)/2) - T 9 (d/2) and T p ((r + d)/2) - T p (d/2) must 
have a common root a = 2 cos ip . Let also set a + r = 2 cos ip . We have 

cospVo = cospi^o and cosgVo = cosg^o- 

The solution is : 

2m7r 2mr 

ipo = ±ip H = T^o H 

P 9 

Since r = 2 cos V'o — 2cos(/?o, excluding the solutions r = which is not acceptable, we 
have the different solutions: 

r = ±4smmir/psmn7r/q, 0<2m<p, < 2n < q . 

It is easy to verify that these results agree with the explicit solutions of the (2,3), (4,3) 
and (4, 5) models. They show that, as expected, all unitary models lead to non-Borel 
summable topological expansions because all terms of the series have the same sign. These 
models thus suffer from the same disease as the pure gravity model. Note finally that 
indeed the number of different values of r is the same as the number of operators in the 
minimal (p, q) conformal model. 

5.4- The p = (2m + l)q ± 1 models 

For m 7^ r is solution of more complicated algebraic equations. In the notation of 
previous subsection we still have: 

, 2nn 

Wo = ±^0 H • 

Q 
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Let us set 

«= 5(^0 + ¥>o), /?= i(V>o -<£o), 
then, making a choice of signs 

TL 7T 

/? = — , r = 4sinasin / 3 = 4sinasin(?27r/(/), 

where 




[(p - 2qrZ)a] sin[(p - 2gZ)/3] = . 



We note that sin[(p — 2ql)(3] = sin(n7rp/g) which can be factorized. We thus find an 
equation for a: 

= sin[(p-2<z0«] = 0. 

This function satisfies the differential equation 

pA(a) (cos qa)' — qA'(a) cos qa = K(p, q) cos a , 

where K is a constant. This equation implies that A(7r/2ij)A(37r/2g) < and thus A(a) 
vanishes at least once in the interval (0, n). We conclude that for all these models the 
topological expansion is not Borel summable. 
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